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Using a new time-dependent measure, we demonstrate for the first time that each defect in a 
representative defect-mediated spatiotemporally chaotic system is associated with one to two degrees 
of dynamical freedom. Furthermore, we show that not all dynamical degrees of freedom are related to 
the defects; additional degrees of freedom are due to underlying phase turbulence. These results yield 
a deeper understanding of the dynamical role of defects and provide hope that these complicated 
systems might be reduced to simpler descriptions. 



05.45.+b, 47.52.+J, 47.54.+r 



I. INTRODUCTION 



The motions and interactions of topological defects 
dominate, at least visually, the dynamics of a wide variety 
of nonequilibrium systems [|lj . This behavior is observed 
in systems as diverse as chemical reactions in shallow lay- 
ers (e.g., §)), thermal convection in horizontal fluid layers 
(e.g., U), aggregation in slime mold colonies (e.g., 
periodically shaken layers of sand (e.g., Jp[), and electri- 
cal activity waves in heart tissue (e.g., Researchers 
have proposed that the macroscopic behavior of these 
complicated systems can be understood by considering a 
simpler system of interacting defects instead of the 
detailed equations describing the system. 

Some nonequilibrium systems containing defects are 
also spatiotemporally chaotic, in that the spatially and 
temporally disordered dynamics continues indefinitely. 
These systems are often termed "defect-mediated" to em- 
phasize the role of the defects. For such systems one 
might intuitively expect that a measure of the average 
complexity D of the system (roughly the minimum num- 
ber of degrees of freedom to describe the dynamics) is 
proportional to the average number of defects {rid{t))t in 
the system. More generally, we consider the possibility of 
a contribution unrelated to the defects, and we account 
for the time-dependence of the quantities: 



D(t) =D b (t) + D d n(t), 



(1) 



where is the average contribution per defect. The 
background term Db(t) is expected to be non-zero if the 
field between the defects does not simply respond pas- 
sively to the motion of the defects. 

Eq. (Q) is important to understand for several reasons. 
First, insight will be gained into what drives the dynam- 
ics of defect- mediated chaotic systems; i.e., whether the 
defects alone determine the dynamics. This result may 
be an important step toward understanding not only the 
building blocks of extensive chaos |j, but also what is 
necessary to control the chaos in these systems. Sec- 
ond, we may begin to understand what determines the 



defect statistics and, in turn, the long-wavelength char- 
acteristics, of these systems. Third, if Eq. (|j) is found 
to be valid, then separate, but coupled, equations for the 
defects and the background might allow a simpler de- 
scription of the complete system Js|,[|]. This description 
might be further simplified through a Langevin equation 
approach in which the fluctuating background field is re- 
placed by a noise term [|l0| — |T"^| - 

We have studied Eq. ^) using a representative defect- 
mediated spatiotemporally chaotic system, the two- 
dimensional complex Ginzburg-Landau equation CD]: 



dtu(x, t) = u + (1 + ici)V 2 u — (1 — iC3 



|m| 2 m, 



(2) 



where u(x, t) is a complex- valued field of size L x L with 
periodic boundary conditions, and c\ and C3 are real- 
valued parameters. Eq. (|^) is an experimentally rele- 
vant amplitude equation universally valid near the on- 
set of a Hopf bifurcation from a stationary homogeneous 
state to an oscillatory state. It is closely related to equa- 
tions describing chemical reaction-diffusion systems and 
excitable media such as heart tissue. For a range of values 
of the parameters (01,03), solutions to Eq. (^J) exhibit a 
defect-mediated spatiotemporally chaotic state in which 
pairs of spirals of opposite topological charge are contin- 
ually created and annihilated Jl3|,p^|. 

For a wide range of values of (ci,c 3 ) in Eq. (g), we 
demonstrate for the first time that the complexity D of a 
defect-mediated spatiotemporally chaotic system is not 
solely due to the defects. In fact, the degrees of free- 
dom can be separated into two sources, suggesting that 
Eq. (]!]) is an appropriate description of defect-mediated 
spatiotemporal chaos. Using a new measure, we show 
that each defect is associated with 1-2 degees of dy- 
namical freedom, and we provide evidence that the re- 
maining degrees of freedom are due to the underlying 
phase dynamics. This result suggests that the chaotic 
dynamics might indeed be modelled using separate, but 
coupled, equations for the defects and the background 
phase-turbulent field. 
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II. MEASURES OF DYNAMICAL BEHAVIOR 

We integrated Eq. (||) using a pseudospectral method 
with time-splitting of the operators, carefully testing for 
convergence with respect to spatial and temporal reso- 
lution and with respect to total integration times. Typ- 
ically, we used a timestep of At — 0.05 time units, a 
system size of L = 64 or L = 32 with 2 Fourier modes 
(or, equivalently, 2 spatial points) per unit of length, and 
total integration times of about 10 5 time units. Defects 
were found by counting the number of crossings of the 
zero contours of the real and imaginary parts of the field 
u{S,t) in Eq. (|). 

The exponential divergence of nearby trajectories in 
chaotic systems suggests that two realizations, u(x,t\) 
and u(x,t 2 ), of a solution to Eq. (||) can be regarded as 
statistically independent when the interval ti — t\ is suf- 
ficiently large. Then, for intervals of duration AT such 
that AT covers many statistically independent realiza- 
tions of the solution, we expect that the root mean square 
deviation of quantities that are each averaged over inter- 
vals of duration AT will converge as AT -1 ' 2 . This con- 
jecture is supported by Fig. |l](a) showing the convergence 
of measurements of the average defect density, n^ AT \t) 
as a function of the duration AT over which each mea- 
surement is averaged. (Each measurement n^ AT ' (t) is 
averaged over the interval t < t' < t + AT.) The mea- 
surements rv- AT \t) converge according to the form 



J AT) 



(<) 



(AT) 



(t) 



1/2 



oc AT 



-1/2 



where () t indicates an ensemble average indexed by the 
time t for the solution u(x,t). 

To characterize the dynamics of the chaotic states, we 
calculated the time-dependent spectrum of finite-time 
Lyapunov exponents (A^ T \t) > A2 (t) > ■ • •) de- 
scribing the average rate of divergence during time inter- 
vals t < t' < t + AT for trajectories in phase space that 
are initially separated by an infinitesimal amount [flB") . 
The finite-time exponents were obtained by an expensive 
but widely-used algorithm [fl6[ that involves integrating 
Eq. (|J) together with many copies of the linearization of 
this equation about the solution u(x,t). 

In previous work on spatiotemporally chaotic systems, 
researchers 0^ |l9| have focused on the infinite-time- 
averaged Lyapunov exponents (and the closely related 
Lyapunov dimension) which average over the fluctuations 
in the dynamics. In this work, we utilize the natural 
fluctuations within the dynamics to extract additional 
information about the system without additional com- 
putational expense. (Because the Lyapunov exponents 
describe behavior in the linear regime about the solu- 
tion, a standard method |16| for obtaining the infinite- 
time exponents is to average the appropriately ordered 
finite-time Lyapunov exponents.) The fluctuations in the 
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FIG. 1. (a) Convergence of the root mean square devia- 
tion of the finite-time defect density rS AT \t)/L 2 as a func- 
tion of the time interval AT. In addition to averages over 
t, results were averaged over 64 different initial conditions for 
(ci , c 3 ) = (3.5, 0.90) in systems of size 32 x 32 with AT = 500. 
(b) Convergence of the root mean square deviation of the fi- 
nite-time dimension D^ AT \t) for a system of size 32 x 32, 
(ci,C3) = (3.5, 1.15), and AT = 32. Solid lines are guides to 
the AT~ 1/2 power law. 



complexity of the dynamics is reflected in changes in the 
spectrum of finite-time Lyapunov exponents. In analogy 
to the Kaplan- Yorke conjecture for (infinite-time) Lya- 
punov exponents we reduce the spectrum to a sin- 
gle quantity by defining a new time-dependent measure 
of the dynamical complexity, the finite-time dimension 



1,(AT) 



(t) 



z?( AT )(<)EEfc( AT )(i)- ]T x^ T \t) 



(\ (AT) (t) 

\ A k(^ T )(t)+l") 



where £;( AT )(i) is the largest integer such that 



v fc <AT) (*) 
Zjj=i 



xi AT) (t)>o. 



Fig. 0(b) suggests that sufficiently separated values 
of this new measure can be regarded as sta- 

tistically independent. We also note that in the limit 
T — > co, D^ AT \t) approaches the (infinite-time) Lya- 
punov dimension D^. All results are obtained in the 
extensive chaos regime for which Dl and {n(t))t grow 
linearly with system size L 2 . 

III. THE DIMENSION PER DEFECT 

Our new measure, the finite-time dimension D^ AT '(t), 
allows the comparison of fluctuations in the complexity 
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to fluctuations in the number of defects (for each set of 
system parameters), effectively separating the degrees of 
freedom correlated with the defects from the remaining 
degrees of freedom. Using our calculation of these time- 
varying quantities, we directly test Eq. ([!]). For a rep- 
resentative set of the system parameters, (c% — 3.5, C3 
1.1.") : . Fig. ||(a) shows the finite-time dimension D^ AT \t) 
plotted as a function of the average number of defects 
n ( AT ) (f ) found during the same time interval of dura- 
tion AT. Each point in Fig. p|(a) represents an average 
of the values of D( AT )(t) for which rv- '(t) is within 
a small range of size 0.2 defects. According to Eq. (Q), 
the slope of the fit yields Dd, the average number of de- 
grees of freedom associated with each defect. The aver- 
age number of degrees of freedom not correlated with the 
defects, (-D&)* in Eq. ([[]), is a significant fraction of the 
the total dimension. (The relatively large contributions 
of {Db)t to the total dimension can provide misleading 
information concerning the dimension per defect if one 
simply compares the (infinite-time) Lyapunov dimension 
to the number of defects.) For a wide range of values of 
finite-time interval AT tested (16 time units < AT < 256 
time units), the values of the slope and intercept are ap- 
proximately constant to within measurement errors, with 
the range of values of Z?( AT '(i) and n^ AT \t) decreasing 
with increasing AT. The accuracy at larger values of 
AT is severely restricted by the decrease in the number 
of measurements and the narrowness of the dimension 
and defect distributions. 

Figures similar to Fig. |2|(a) are obtained for a wide 
range of parameters (01,03). Fig. |2|(b) shows the values 
of Dd obtained from linear fits of the data for each set of 
system parameters. We find that Dd decreases slowly as 
the number of defects increases. We note that a rough 
extrapolation of the data in Fig. ||(b) toward zero defects 
suggests an approximate value of Dd ~ 2. 

Further work (at great computational expense) will be 
necessary to discern the physical origin of the measured 
values of Dd- The value of Dd may indicate the aver- 
age fractal dimension of the trajectories of defects; for 
example, a value of Dd = 2 would be expected if the 
trajectory of a defect is essentially random in two spa- 
tial dimensions. Possible explanations for the decrease 
in Dd at larger defect densities are correlated behavior 
between defects and constraint of the defect trajectories 
due to the proximity of other defects. Evidence of these 
possibilities might be found in the small amount of neg- 
ative curvature observed in Fig. |2|(a) and similar plots 
for other parameter values. However, much more data 
will be needed to measure the curvatures accurately and 
to assess the impact of these small curvatures on other 
quantities reported here. The small magnitudes of the 
curvatures are consistent with analytical and numerical 
work showing that defects in Eq. (||) are weakly interact- 
ing §■ 

We are currently investigating these ideas, and we are 




FIG. 2. (a) Finite-time dimension D i -^ T \t) vs. the av- 
erage number of defects n AT (t) during the same inter- 
val AT = 32 time units for (01,03) — (3.5, 1.15) in a system 
of size 32 x 32. Each point is an average over all measure- 
ments with n AT (t) within bins of width 0.2 defects. The 
solid line is a linear fit with each point weighted by the 
number of measurements in the respective bin. (b) Val- 
ues of Dd for various parameters (01,03). Circles indicate 
ci = 3.5 and 0.85 < 03 < 1.55; triangles indicate 03 = 0.9 and 
2.5 < ci < 5.0. 



studying whether similar values of Dd are found in other 
defect-mediated spatiotemporally chaotic systems. 



IV. THE PHASE CONTRIBUTION 

The non-defect degrees of freedom (-Dfe)t vary signif- 
icantly over our range of system parameters, as seen in 
Fig. ||. We speculate that this "background" contribution 
is due to an underlying phase turbulence. To elucidate 
this further, we obtained an estimate of (Db)t using a 
phase equation derived in the limit v — (c\Cz — 1) — > 
|p0f . Upon rescaling, the phase equation can be written 
in a parameter-less form known as the two-dimensional 
Kuramoto-Sivashinsky equation: 

d t 8(x, t) = -V 2 8 - (V6») 2 - V 2 V 2 6>. (3) 

The scaling leading to Eq. (|J) introduces a factor 
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allowing the computation of a "phase dimension" Dg 
of Eq. (Q) as a function of (01,03) and the system 
area Acgl = L x L, given a single measurement of the 
(infinite-time) dimension Z?ks °f Eq. ([|) for a system of 
size A KS p): 
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FIG. 3. Background dimension (Db)t (triangles) and phase 
dimensions Do (circles) for (a) C3 = 0.9 and 2.5 < ci < 5.0, 
and (b) ci = 3.5 and 0.85 < c 3 < 1.55. 
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Remarkably, as seen in Fig. |^, the behavior of the phase 
dimension Dg as a function of the average number of de- 
fects (n(t))t is similar to the behavior of the measured 
average background dimension (-D&)* (although the pre- 
diction of Eq. (^) is smaller than the measured values). 
This is particularly evident in the non-trivial behavior on 
the right-hand side of Fig. ||(a). This result is particu- 
larly surprising since Eq. (|J) is derived from an expansion 
in the small parameter v = (ciCg, — 1) — > 0, while, for the 
solutions studied here, v — 0(1). Additionally, since 
the defects are not bound in pairs of opposite topological 
charge, the phase of the complex field should be globally 
influenced by the presence of the defects. The effect of 
the defects on the phase dimension might not be as large 
as one would naively expect since chaotic lengthscales 
have been shown to be much shorter than lengthscales 
corresponding to the phase [Q,|l^,[l9); however, this ef- 
fect might account for the increased difference between 
Dg and (Db)t when the system contains more defects, as 
seen in Fig. |^(b). That similarities remain even in the 
presence of these effects suggests that the background 
dimension is closely related to the phase turbulence. 

Additional insight might also be gained by studying 
the solutions to Eq. (^) and the solutions to Eq. (||) for 
some small values of v. These solutions are purely phase 
turbulent cellular patterns |l4| with the cells delineated 
by shocks in the phase field. A cellular pattern also can 
be found in the defect-turbulent state of Eq. (|^) with the 
defects moving predominantly within the cells. Further 
work will be required to test whether the number of these 



cells is related to the background dimension. 



V. CONCLUSIONS 

Prompted by measurements of the ergodic nature of 
spatiotemporally chaotic attractors, we have defined a 
new time-dependent measure, the finite-time dimension. 
By employing this measure we have found a linear rela- 
tionship between the degrees of dynamical freedom and 
the number of defects in the system during the same 
time interval. These results show that each defect in 
a defect-mediated spatiotemporal chaotic system is as- 
sociated with Dd w 1-2 degrees of freedom. However, 
a large number of degrees of freedom are not associated 
with the defects, and we have provided evidence that 
these degrees of freedom are related to the underlying 
phase turbulence. Our results yield a deeper understand- 
ing of defect-mediated chaos and provide hope that these 
complicated systems might be reduced to simpler descrip- 
tions. 

We thank E. Bodenschatz, I. Melnikov, and H. Riecke 
for valuable discussions. This work was supported by the 
US National Science Foundation and the Cornell Theory 
Center. 
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